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Abstract 

We propose a novel cosmological scenario, in which standard inflation is replaced by an expanding 
phase with a drastic violation of the Null Energy Condition (NEC): H 3> H 2 . The model is 
based on the recently introduced Galileon theories, which allow NEC violating solutions without 
instabilities. The unperturbed solution describes a Universe that is asymptotically Minkowski in 
the past, expands with increasing energy density until it exits the regime of validity of the effective 
field theory and reheats. This solution is a dynamical attr actor and the Universe is driven to it, 
even if it is initially contracting. The study of perturbations of the Galileon field reveals some 
P^ | subtleties, related to the gross violation of the NEC and it shows that adiabatic perturbations are 

cosmologically irrelevant. The model, however, suggests a new way to produce a scale invariant 
spectrum of isocurvature perturbations, which can later be converted to adiabatic: the Galileon is 
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forced by symmetry to couple to the other fields as a dilaton; the effective metric it yields on the 
NEC violating solution is that of de Sitter space, so that all light scalars will automatically acquire 
a nearly scale-invariant spectrum of perturbations. 



1 Introduction 
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The connection between early cosmology and high energy dates back to the discovery of the ex- 
pansion of the Universe. Going backwards in time, the Universe contracts and the energy density 
increases more and more until it eventually reaches Planckian values and General Relativity breaks 
down. The underlying assumption to this general argument is that the stress-energy tensor T^ v 
satisfies the NEC, which states that T iw k li k v > for every null vector k^. For a perfect fluid this 
is equivalent to the inequality p + p > 0, where p and p are respectively the energy density and 
the pressure. For a Friedmann-Robertson- Walker (FRW) metric this inequality implies that the 
energy density, and therefore the Hubble parameter H, decreases as the Universe expands, as the 
covariant conservation of the stress-energy tensor reads p = —3H(p + p). 

If a violation of the NEC were possible, then a Pandora's box of non-standard cosmologies 
would open up and in particular the contraction of the Universe going backwards in time would not 
necessarily lead to higher and higher energy densities. The realization that energy conditions can 
be violated, although "standard" matter satisfies them, has a notable history in cosmology. Indeed 



the strong energy condition, equivalent for a perfect fluid to the inequality p + 3p > 0, implies 
that the expansion of the Universe is always decelerating a < 0, which resonates well with the 
Newtonian intuition. The evidence of the present acceleration, however, strongly indicates that the 
fluid which now dominates the Universe violates the strong energy condition and the same happens 
in the past during inflation, the most compelling theory of the early Universe. Given that these 
two important revolutions in cosmology are based on a violation of the strong energy condition, it 
is natural to wonder whether we are missing something taking the NEC for granted. It is worth 
emphasizing that the NEC is usually taken for granted not due to lack of imagination, but because 
of its exceptional robustness — it is especially hard to construct consistent effective field theories 
that violate it pQ. Morevover, perhaps less decisively, the NEC protects standard properties of 
black-hole theormodynamics. 

In this paper we describe a novel cosmological scenario in which the NEC is grossly violated, 
and the Universe starts from a very low energy state, asymptotic to Minkowski in the far past. 
As pointed out in [2], the violation of the NEC is possible in the context of the recently studied 
Galileon theories [3]. For these theories the usual relation between the violation of the NEC and 
the presence of pathological instabilities [lj is avoided, due to the presence of higher derivative 
interactions. A similar situation happens in the context of Ghost Condensate theories [HE], with 
the important difference that here we will have a strong violation of the NEC, i.e. H S> H 2 , while 
in the previous models, only a moderate violation of the NEC is compatible with the stability of 
the system [5]. 

A general class of Galileon theories, as we will see in Section [21 gives rise to a very peculiar 
evolution of the scale factor with the Hubble parameter becoming larger and larger as the Universe 
expands: H oc (— £)~ 3 as t — > 0~. This implies that most of the energy is created suddenly, with 
the scale factor blowing up as a ~ exp(l/t 2 ), in a sort of Genesis which (partially) justifies our 
dubbing of the scenario. As p increases, the system will eventually exit the regime of validity of 
the Galileon effective field theory and here we assume that the energy gets transferred to more 
conventional degrees of freedom in a reheating process, similarly to what happens in inflation. 
This background evolution is completely stable and it represents a dynamical attr actor. Notably 
the Universe evolves to this expanding phase even if it is initially contracting, a behaviour which 
is only possible because of the violation of the NEC. It is remarkable that this scenario in some 
sense explains why the Universe is now expanding, while we are usually forced to postulate initial 
conditions with a large positive H , which then goes on decaying for the entire evolutional 

The study of perturbations in Section [3] and Appendix [A] shows various peculiarities of this 
model. One such peculiarity is that the energy density of the background solution vanishes in 
the limit in which gravity is decoupled, p oc l/M™. Another peculiar feature is that the leading 
adiabatic solution does not correspond to the standard £ = const, mode, but to a constant time 
shift of the unperturbed solution (this will be shown in Newtonian gauge in section U]). Such a 
mode is going to decay during the standard post-reheating FRW evolution and this will allow us 
to conclude that the fluctuations of the Galileon field do not give rise to any relevant cosmological 
perturbations on large scales. Actually in Appendix[B]we will show that the Galileon perturbations 
are not amenable to any classical interpretation as they experience no relevant squeezing. 

Fortunately, another source of scale invariant perturbations is naturally present in our model, 



1 Oi course our 'auto-expanding' solution run backwards in time is also a solution — i.e. we also have 
solutions that contract for all times and approach Minkowski space in the future. What we want to stress 
is that we are driven towards our expanding solution starting from an unusually large basin of attraction, 
which includes contracting initial conditions as well. The time-reversed solutions we just alluded to start 
outside this basin. 



as we explain in Section^ The Galileon Lagrangian is invariant under the conformal group SO(4,2) 
and the time dependent solution breaks it down to SO(4,l), the isometry group of de Sitter space. 
The only way another field can couple to the Galileon while respecting the conformal symmetry is by 
treating the Galileon as a dilaton, that is through a fictitious, conformally flat, metric. Therefore 
all other fields will perceive the Galileon background as a "fictitious" de Sitter space and their 
dynamics will be essentially the same as for inflation, even though the Einstein metric at the time 
when cosmological perturbations are generated is virtually flat. In particular a massless scalar will 
acquire a scale-invariant spectrum of perturbations. These isocurvature perturbations can then be 
converted to adiabatic in a variety of ways, exactly as it happens for inflation. This novel mechanism 
to produce a scale invariant spectrum of perturbations shares some similarities with the attempts 
to explain the present acceleration through the universal coupling of matter to a scalar field. In 
both cases, an approximate de Sitter space is realized not in the Einstein metric but in the Jordan 
one. However, it is crucial to keep in mind that this "fake", Jordan-frame de Sitter space seen by 
fluctuations, is by no means helping us solve the horizon problem and that the peculiar cosmological 
history we outlined here and that we are going to describe at length in the paper happens in the 
Einstein frame. Indeed our system violates the NEC even in the absence of dynamical gravity 
[2]; if we do turn on dynamical gravity, coupling it minimally to our system, this NEC- violating 
stress-energy tensor will generate an Einstein- frame NEC-violating geometry. Related to this, we 
also notice that the Galileon was originally introduced as a possible explanation of the present 
acceleration [3] , benefiting from its natural implementation of the Vainshtein screening mechanism 
at short scales. Here we are using it in a completely different context, motivated by its healthy 
violation of the NEC. 

It is important to stress a problem with our model: superluminality. Perturbations around the 
SO(4,l) invariant background move at the speed of light due the large amount of residual symmetry, 
and actually gravity corrections to the solution make them slightly subluminal. On the other hand, 
if we allow for large departure from the background, perturbations around the new solution will 
move superluminally. In Section [6l we will discuss this issue and its implications applying to our 
case the general discussion of [2] . Conclusions are drawn in Section [7J 

2 The background 

Our starting point is the simplest version of the conformal Galileon minimally coupled to gravity: 
the Lagrangian for the scalar field is just the sum of the kinetic term and the Galilean invariant 
cubic interaction plus the (dir) term needed to recover conformal invariance |3j 
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Lorentz indices are contracted with g^ v and the □ contains a covariant derivative^. Notice that 
the conformal symmetry of the ir Lagrangian is explicitly broken by the coupling with gravity. 
We could add all Galilean-invariant interactions together with their conformal completions [3] , and 
in fact a fully consistent NEC-violating system obeying all requirements of [2] will have them. 
However our analysis and results below would not be affected in an essential way, since virtually 
all our results follow from the symmetry structure of the theory. One important difference is that, 
with our minimal Lagrangian, the kinetic term has the wrong, ghost-like sign, around the trivial 
background it = 0, while it is healthy in more general Galilean Lagrangians [2J. This instability 



2 We are using the mostly plus signature. 



is not relevant for us as we are going to be interested in a different background solution, but it 
may become important if the system eventually evolves to tt = after reheating. In the present 
analysis we stick to the minimal theory ([1]), for the simplicity of the computations involved. We 
will comment on the effect of adding higher order Galilean terms when relevant. 

The signs have been chosen so that if gravity is decoupled this action has a solution in Minkowski 
spacetime of the form 

e^s = __L _oo < t < o , (2) 

provided that 

H 2 - 2A3 ™ 

H ° - w • (3) 

Such a solution spontaneously breaks the conformal group 50(4,2) down to the de Sitter group 
SO (4, 1). More importantly for our purposes, this "de Sitter" field configuration violates the NEC 
but has no instabilities [2]. The tt stress-energy tensor can be easily computed from the action (TjQ) 

ca T — 2 (55V 

aS ^ - y/=gS 9llv > 

V = -/ 2 e^[2cVcV-V,(37r) 2 ] 

f 3 
- ^ [2 d^d.irOiT - (cV d^Orr) 2 + d„ir d^dirf) + 9{a , d a ir d a (dirf] 

'" [Aidnfd^d^-g^idTT) 4 ] . (4) 
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By plugging the solution ([2]) into this expression with g^ u = r]^ u we find that it has vanishing 
energy density — this is a consequence of scale-invariance which is left unbroken by the background 
[2] — and negative pressure oc — 1/t 4 . 

We defer a thorough stability analysis for this solution until the next section. For the moment, 
let us consider the dynamics of homogeneous perturbations 5n(t). From the Lagrangian ([1]) we 
immediately get the equation for 5n in the linear regime: 

2 4 

5if 5tt ^5n = ; (5) 

the two independent solutions are 5tt ~ 1/t and 5it ~ t 4 . The latter decays away for t — > 0~ and is 
thus no source of worry. The former blows up at late times, but in fact it just describes the same 
background solution we are interested in, slightly shifted in time [2]: 

TidsO + e ) - 7Tds(i) + 7i"dsO) • e = 7Tds(i) - - • (6) 

We conclude that a generic homogeneous initial condition that corresponds to a small departure 
from 7Tds will be diluted away at late times. 

Now let's reintroduce the coupling to gravity: the presence of a non-zero stress-energy tensor 
will source a gravitational field and tt = tt^s with flat metric will no longer be a solution. However, 
since also the pressure vanishes at early times (t —> — oo) the gravity- free solution is recovered in 
this limit, while as time goes on corrections to this asymptotic behavior will become larger and 
larger. As we are interested in cosmological solutions of the form 

ds 2 = -dt 2 + a 2 {t)dx 2 , 7r = 7r (t) (7) 



we have just to solve Friedmann's equations for the Hubble rate H , 
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and the it e.o.m. will be automatically satisfied. 

From eq. (JD) we get the energy density and the pressure, 
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(11) 



where we used eq. ([3]) to write A 3 in terms of Hq. We cannot solve Friedmann's equations ana- 
lytically. Still, we can compute the asymptotic behaviors of the solution at early and late times. 
As we said, at early times the solution is approximately the unperturbed one tv = ir^s and H = 0, 
with small corrections proportional to G = 2 , which we can calculate perturbatively. Since 

p ~ 0(G) <C p = Pds + C(C?), the leading contribution to eq. ([9]) is H ~ — AirGp^s and this 
expression can be integrated to find the Hubble rate at early times: 
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This result has a number of unusual properties: i) we can add to the value of H an arbitrary 
integration constant and still have a solution of eq. ([9]); we will discuss this possibility shortly, 
while for the moment we set the constant to zero; ii) because p — > in the limit Mp\ —> oo, the 
Hubble rate is proportional to 1/Mpj, unlike for standard cosmological scenarios where it scales like 
1/Mpi; Hi) H increases with time as a consequence of the NEC violation, with a rate H 3> H 2 . 

Having computed the value of H we can plug it into ([8j) , or equivalently into the scalar equation 
of motion, and extract the 0(G) correction to tt^s 
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(13) 



(we are choosing ttq — > ir^s for t — > — oo as initial condition, as required by consistency of the 
approximations we have adopted so far.) At late times, for t 2 < -772-777 the above approximation 
breaks down. Numerically integrating the Friedmann's equations shows that both ir and H diverge 
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their asymptotic behaviors: 
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which indeed match the actual numerical solutions. This gives the peculiar evolution of the scale 
factor 



a(t) ~ exp 



1 



3H 2 M 2 { (t - t) s 



(16) 
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Figure 1: The cosmological evolution in our model. 

We have then the following scenario, represented in figure [TJ The Universe starts at t — > — oo 
in a quiescent state, with flat metric and ir = vr t jg(t). Asymptotically in the past this configuration 
has zero stress-energy tensor, and so it is a solution. Then, as time goes by, a negative pressure 
arises, p ~ — 1/t 4 , which makes the Universe start expanding and the energy density grow. The 
actual solution then departs from the original "de Sitter" configuration. H, p, and p grow more and 
more, until 7r becomes strongly coupled. At that moment the effective theory of it breaks down, 
and we cannot predict what happens next. We can imagine that at that point, most of the available 
energy gets converted into radiation, the Universe reheats, and the standard radiation-dominated 
era takes over. In this era the Galileon 7r may evolve to the n = const, solution, or cease to be a 
good degree of freedom. 

The strong-coupling scale of the theory "runs" with = Hc ) e 7r , and it is [2] 
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9 = H /f . 



(17) 



Notice that this estimate does not take into account gravity and the associated explicit breaking of 
conformal symmetry. It may thus get modified when dynamical gravity is included. Our effective 
theory breaks down — at the latest — when typical energies become larger than A str0 ng- A good 
measure of a cosmological solution's typical energy is the freeze-out frequency for fluctuations: if 
fluctuations are strongly coupled at freeze-out, the background solution is hardly consistent. In 
the phase where things are blowing up, eqs. (fbfl[T5]) . freeze-out happens at frequencies of order H. 
The highest Hubble rate we can get before strong-coupling/reheating thus comes from equating 
eq. (fT5|) with the strong-coupling scale (fTTj) . We get an impressive 
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(18) 



at which GR breaks down anyway. Because of this, we don't want to go this high in H, so we will 
assume that our effective theory breaks down before becoming strongly coupled, which is also very 
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much in line with the arguments presented in [2j. Notice that all other measures of the solution's 
typical energy — for instance 4>/4> etc. — are smaller that the one we have adopted, and are thus less 
constraining. In conclusion, the reheating temperature is essentially unconstrained in our model. 

At this point one may ask whether the background discussed so far is an attractor. As in 
the previous discussion in the absence of gravity, we can study the homogenous perturbations 
7r(i) = 7ro(t) + 5ir(t) = 7Tds(i) — \4tt (n't) 2 ~*~ 'W*)- To expand linearly in 5tt the only conditions 
needed are 5ir <C ttq and Sir <C 1. However even in the linearized approximation there are two 
possible regimes, depending on whether the perturbations are also smaller than the gravitational 
corrections suppressed by l/Mp^ or not. 

Let's start from the scalar equation of motion: 



e 27r 



n + -k 2 ) - -^n 2 n = -3e 27r H*fr + \ [AH-frJf + 2ir 2 (3H 2 + H) + 2Hir 3 ] . (19) 



Keeping 5ir fixed and going to early times we can neglect gravity corrections and the linearized 
equation of motion for Sir reduces to the one studied before, eq. ©, with the two solutions {1/t; t }. 
Notice another unfamiliar feature of this background: since pas = and we are in a regime where 
gravitational contributions are small, linear perturbations to the scalar give a contribution to H 

larger than that of the background. 

f 2 i 
In the opposite regime 5tt <C -fW- , H t \i terms proportional to H, H on the right-hand side of 

(|19p will give a contribution to the linear equation for the pertubations. Since H 2 , contrary to 



the previous case, is now dominated by the background solution we can expand in equation (jHJ) 
5H 2 = 2H5H together with 7r = ttq + 5ir to find SH = —5ir/t — 5-k. Using this expression in the 
RHS of (fT9|) gives the perturbations' equation in this regime: 



5tt Sir 

T~T 2 



^ + — - 72" = ° • ( 2 °) 



The two solutions are {1/t; £}; we have again the shift in time and a decaying solution that implies 
convergence to the attractor. 

We can now study the most general solution for the Hubble rate. Suppose we start with an initial 
perturbation Sir = At 4 "; this gives a constant contribution to the energy density p = —10Af 2 /HQ, 
which corresponds to the integration constant for H we alluded to below eq. (|12p . Since this can be 
positive or negative, a generic initial condition can produce a background that starts in expansion 
or even in contraction. 

No matter which sign we choose for the initial condition, the perturbation decays as time goes 
on and we are eventually driven to the original expanding background. Notice that as we approach 
the unperturbed solution, 5ir will move from the first regime — where it dominates the energy 
density and has 5ir ex t 4 as solution — to the second one, where 5ir is smaller than the perturbation 
induced by gravity. In the intermediate regime we have no analytic control and one may be worried 
that we do not recover the unperturbed solution eventually. However we know that the equation 
H ~ — 47rGpdS is always a good approximation since p <C p- It tells us that H follows the evolution 
discussed above up to small corrections, even in the transition between the two limiting regimes 
for Sir, where we don't have an explicit solution for homogenous perturbations. 

In conclusion the NEC violating solution is an attractor for general homogeneous initial condi- 
tions close to the de Sitter solution: 5n <C vr^s- 



3 Scalar perturbations 



Let us now move away from homogeneity and discuss scalar perturbations. To begin with, let 
us first assume again that gravity is decoupled, Mpi — > oo. Since the solution ([2j) spontaneously 
breaks the conformal group SO(4, 2) down to the de Sitter one 5"0(4, 1), the Lagrangian for small 
perturbations will be invariant under the de Sitter symmetries, whereas the broken symmetries will 
be non-linearly realized. In particular the fluctuation £(x) defined via n(x) = 7Tds(i + £,(%)) is the 
Goldstone boson associated with the spontaneously broken time-translational invariance t — >• t + e, 
which is now realized non-linearly as £ — > ^ + e. Indeed from eq. ([T]) we get the quadratic Lagrangian 
for £ 






^ = -^-MY , (2i) 



which is manifestly shift-invariant. The kinetic energy is positive, thus ensuring stability for the 
background solution against short-wavelength perturbations. For long- wavelength ones instead we 
get back to eq. ©, now written in terms of £ : 

L - ~ t ik = , k < l/t . (22) 

The solutions are 

£ fc ~ t 5 , const . (23) 

The constant one dominates at late times and simply describes — now manifestly — the original 
background solution slightly translated in time. We thus conclude that, in the absence of gravity, 
the solution 7Tds is an attractor also for initial perturbations with non-vanishing gradients. 

We now turn on gravity and in general we expect that the dynamics of scalar perturbations at 
large distances will be modified by their mixing with the scalar sector of g^ u . Let's see how this 
works explicitly. Suppose we have the background solution wo(t), H(t). Then, if we consider small 
fluctuations, it is particularly convenient to work in 'unitary gauge': 

7r(x,t) = 7r (t). (24) 

This fixes time-diff invariance: we are defining the equal-time surfaces as the equal- n ones, and the 
pace of time as that of the unperturbed solution. In this case there are no fluctuations in n, and 
the scalar fluctuation is in the metric tensor. We will fix the space diffs later. 

Following [6] and [5], we use ADM variables for the metric: the induced 3D metric gij, the lapse 
N = 1/a/— 5 00 , and the shift N{ = g^. It is straightforward (see Appendix |A|) to write the full 
action, the Einstein-Hilbert one plus the Galileon part, using these variables: 

S g = ^M^J^x^NlRz + iKijKV-K'i 2 )} (25) 



S* = f 2 Jd 4 x^N 



6 ° iV 2 + 9if 2 iV 3 l + m 2 N\ 



(26) 



where and henceforth spatial indices are raised and lowered via the spatial metric gij , and we have 
used the extrinsic curvature of costant-t hypersurfaces 

K * = Jn \- dt9i3 ~ ViNj ~ VjNi \ ■ (27) 



In fact the structure of the action is largely constrained by symmetry considerations [5j. The 
background solution spontaneously breaks time translations (cliffs) as well as Lorentz boosts. This 
is made explicit by working with ADM variables, and, in unitary gauge, by allowing for explicit 
functions of time in the action. Then, we just have to write down all possible operators compatible 
with the residual symmetries, namely time- and space-dependent spatial diffs, x % — > x l + £ % (t,x), 
and 3D rotations. The generic Lagrangian for matter (w, in our case) then is [5] 



£r = Jd 4 x^N 



+ -M 4 (t)(5N) 2 -M 3 (t)^i(5iV + 



(28) 



The terms in the first line are the only 'tadpoles' there are: they start linear in the metric pertur- 
bations, thus yielding a non-trivial stress energy tensor on the background solution. As a result 
their coefficients are uniquely determined in terms of the background H (t) by the Friedmann equa- 
tions. The terms in the second line start quadratic in the fluctuations, and their coefficients are 
unconstrained. 6N is obviously the fluctuation in N, N = 1 + 5N, whereas the tensor SEij is, apart 
from an extra factor of N, the fluctuation in the extrinsic curvature of constant-t surfaces, 

Eij = NK tj , SE id = Eij - a 2 H 9ij . (29) 

Finally, the dots stand for higher-derivative terms — which in our case vanish, because of the magic 
properties of our conformally invariant Lagrangian — and for interaction terms, cubic and higher in 
the metric fluctuations — which we are not interested in. At the quadratic, two-derivative level the 
action (f28|) is all we need. 

Our action (126p can indeed be recast in the form (128p (see Appendix |A]) with 

M4 W = 14 (24 + 4*0 + 9#4) , M 3 (t) = |^4 . (30) 

In the following, however, we keep the analysis as general as possible, because then we can apply 
it to other conformally invariant Lagrangians as well. However from the expression for M 3 above 
we see why we can violate the 'theorem' of ref. [5]: there it was assumed that the rate at which the 
Lagrangian coefficients — in particular M 3 (i) — vary with time is at most of order H . Here however 
at early times H ~ ^/(M^H^t 3 ), whereas (l/M 3 )<9 t M 3 ~ 1/i, which is much larger than H. 
Thus our example does not satisfy the hypotheses of the theorem. On the other hand, at late 
times the rate of M 3 {t) is slower than the Hubble rate. However at late times H <C H 2 , which, 
according to ref. [5] is compatible with a ghost-free violation of the NEC. In conclusion: there is 
no contradiction with our NEC-violating system being free of instabilities throughout. 

Before proceeding, it is time to comment on what changes if we allow for a more general 
conformal Galilean Lagrangian to start with. The NEC violating background solution will be the 
same [2] , apart from numerical factors and this fixes the first line of (|28|) . Galilean theories give rise 
to equation of motion containing at most two derivatives on each field |3j , so that at the quadratic 
level the operators (5N) 2 and SE^SN are the only possible ones: all the others would give rise to 
equation of motion with more than two derivatives on £ [5] . With a proper choice of the coefficients 
of the Galilean Lagrangian the fluctuations around the NEC violating solution are healthy, as in 
our example eq. ([1]), and it is also possible, at the same time, to have stable perturbations around 
the 7r = background [2], i.e. to flip the worrisome sign of the first term of eq. (Q]). Also the time 
dependence of eq. (|30p will remain the same in the first phase |t| 3> to f° r a general Galileon theory, 
so that all our conclusions can be straightforwardly applied to the more general case as well. 
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We finally move to compute the quadratic action for the propagating scalar mode. This was 
already done in ref. [5] for the Lagrangian (|28p . but only in the M ,M 3 = const case, which as 
we argued is quite different from ours. Following Maldacena [6], the spatial cliffs can be fixed for 
instance by imposing 

9ij = a 2 (t) [(1 + 2C)6ij + 7ii ] , d, nij = , ju = . (31) 

The transverse traceless matrix 7y corresponds to tensor modes, which we will discuss below. 
For the moment we can consistently set jij = 0, since 3D rotations are left unbroken by the 
background solution, thus preventing any mixing between scalar and tensor modes at the quadratic 
level, 2(8)0 75 0. £ parametrizes the only scalar propagating d.o.f. As we will now see, the remaining 
metric components, 500 and goj, can be expressed as functions of £ through the constraint equations. 
These are the variations of the full action S = S g + S n with respect to N and N° (see their explicit 
form in the Appendix [A]). At zeroth-order in the fluctuations, the constraints are solved by the 
background solution. In particular, the Hamiltonian constraint reduces to Friedmann equation, 
whereas the momentum constraint is trivial. To get the quadratic Lagrangian for £, we need to 
solve the constraints at first order in the perturbations £, 5N, and N J . Defining N J = dj/3 (we can 
set to zero the transverse vector piece in N J , for the same reason as for the tensor modes), we get 

SN = iMgM ( 32 ) 

V72/3 _ 2M l\ 1 nV 1 -4Afg 1 i/-12Mg 1 HAf 3 +3Af 6 +2Mg 1 M 4 l , ■. 



The quadratic action for £ is then obtained by plugging these back into the original action, eq. ([25 
plus eq. (J28j) . After some integrations by parts we get 



S f 



Jd 4 xa 3 [A(t)C 2 -B(t)^{V() 2 } , (34) 



where 



Ml ( - 4M 4 H - 12Af2, HM 3 + 3M 6 + 2Ml M 4 ) , , 

A(t) = — ^ : — * (35) 

(2M 2 l H-M^) 2 

Ml(- 4M 4 H + 2Ml HM 3 - M 6 + 2Af2, d t M 3 ) , s 

Bit) = — ^ ^ ^ — = ^— '- . (36) 

(2M^H-M3) 2 

As a check, notice that for a cosmology driven by a minimally coupled scalar with a standard 
kinetic term and non-derivative interactions, we have M 3 = M 4 = 0. This implies A(t) = B(t) = 
-M^(H/H 2 ), which is Maldacena's result [§]. 

At early times — or equivalently at leading order in 1/Mpj — we have (see eqs. (|12|) and (|30p ) 

H~--J- — m 4 ~ — — M 3 i^_l l+l^iff-i (17) 

so that A(t) and !?(£) above reduce to 

A(t) = B(t) = 9 -^lt 2 . (38) 

J 2 
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Also, given the smallness of H, the scale factor can be approximated as constant, a{t) = l + 0(l/i 2 ). 
Therefore at early times we have 

5 C = ^/^(^) 2 [C 2 -(VC) 2 ] • (39) 

It is easy to deduce the spectrum of Q directly from this action. The action is invariant under 

1 
t,x-> At, Xx C -> To C ; (40) 

A^ 

as a result, the equal-time 2-point function of £ must have the general form 

(C(t,0)C(t,x)) = -^^F(\S\/t), (41) 

where F is a generic function, with no additional dependence on the model parameters (the nor- 
malization prefactor comes from the overall constant appearing in the action (|39p .) Notice that it is 
crucial that we look at symmetries of the action rather than simply at symmetries of the equation of 
motion (under which the action might change by an overall multiplicative constant). Indeed, if we 
want the n-point function to have the same transformation properties as Q n , the vacuum state has 
to be invariant under the symmetries considered. This is the case, barring spontaneous breaking, 
if the action is invariant. 

At short distances, \x\ <i, we have to recover the standard Minkowski 2-point function for a 
(non-canonically normalized) massless field: 

\x\ 2 
F (\x\/t) ~ -£- , W\ « 1*1 • (42) 

To get the behavior of F at large distances, we use the fact that the quantum ( solves the classical 
equations of motion. For long wavelengths, or equivalently late times, we have the two behaviours 

C ~ const , - . (43) 

The second term dominates so that (CC) ~ l/£ 2 also at late times. This implieqj 



3 Using the same method for inflation, one would start with the action in conformal time rj of the form [6] 

S = M> { e J d*x j^ [C' 2 - (Vtf] , (44) 

where H and the slow-roll parameter e can be taken as constants at leading order in slow-roll. This action 
is invariant under 

r],x —> \r],\x C ->• ( , (45) 

which implies 

(C(r ? ,0)C(r ? ,f)) = ^F(|^/r 7 ), (46) 

with F ~ rj 2 /\x\ 2 for |x| <C \rj\ to reproduce the Minkowski result and F ~ const for \x\ ^ |?y| to reproduce 
the time evolution £ ~ const one deduces from the equation of motion. In Fourier space this gives the 
celebrated 1/k 3 spectrum. 

11 



which is a very blue spectrum, going as k 1 in Fourier space. Indeed the standard calculation in 
terms of modes (see Appendix [B]) gives 

<«'•*>«*.»» = <*0''< S+ *>]IiSl*5P ' <48> 

Although the spectrum clearly shows that the Galileon perturbations are irrelevant on large 

scales, the reader may be puzzled. Indeed, the two independent solutions of the equation of motion 

are 

sin kt cos kt 

~kT ' ~kT { > 

which, in the limit of long wavelength, respectively give as we discussed 

£ ~ const , - . (50) 

The first solution describes the celebrated conservation of £ on super-horizon scale^l However, this 
solution is irrelevant as the second mode dominates at late times. On the other hand we saw that, 
in the absence of gravity, we have an attractor and we thus expect to flow to the adiabatic solution 
C = const. To understand why it is not the £ = const mode that dominates, in the next Section 
we move to Newtonian gauge where things will (hopefully) clarify. The Section may however be 
skipped without loss of continuity. 

As to the cosmological (ir)relevance of adiabatic perturbations, in Appendix [B] we study the 
quantum state of each Fourier mode until it eventually comes back into the Hubble radius, cal- 
culating the amount of squeezing induced by the cosmological evolution and comparing our case 
with inflation. We will see that no appreciable squeezing is produced and all Fourier modes are 
practically in the ground state when relevant for observations. This unambiguously shows that 
there are no sizable perturbations. 

4 The two adiabatic modes in Newtonian gauge 

To clarify the situation it is better to write things in a way that has a smooth limit when gravity 
is decoupled, which is clearly not the case in the gauge we are using, where perturbations of the 
scalar are set to zero. Things are much clearer in Newtonian gauge 

ds 2 = -(1 + 2$>)dt 2 + a 2 {t){\ - 2^)dx 2 . (51) 

Instead of doing an explicit change of gauge we can equivalently use the Bardeen potentials (for 
a recent review on cosmological perturbation theory see [9]), i.e. the gauge invariant combinations 
which coincide with the scalar perturbations in Newtonian gauge. These are expressed in terms of 
our variables as 

$ = 5N+(a 2 p) (52) 

tf = -C-Ha 2 f5. (53) 

Also the perturbation of the scalar field in Newtonian gauge can be written as a gauge invariant 
combination which in our notation reduces to 

£ = a 2 /3 . (54) 



4 See [8] for a proof of the conservation of £ at all orders in perturbation theory, that applies to the 
Lagrangian studied here. 
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First of all, let us verify that in this gauge one has a smooth Mp\ — > oo limit |j. We have to 
check that in this limit, keeping the amplitude of the scalar mode perturbation £ fixed, the metric 
becomes Minkowski, i.e. <J> and \& go to zero. From eq. (|33|) one can see that in the Mp\ — > oo limit 
with fixed £, which is equivalent (via eq. ([51]) ) to fixed /?, £ goes to zero. As H also goes to zero, 
ty in eq. (|53p goes to zero. For <3? the limit is not so evident as both terms in eq. (|52p remain finite. 
However one can check that they cancel in the limit of decoupled gravity, by using the equation of 
motion of £ derived from the action (|39p. As spacetime approaches Minkowski when Mp\ — > oo, 
also the equation of motion for the scalar will reduce to that in the absence of gravity. 

We are now in a position to compare with the homogeneous perturbations studied in Section 
[2j Rewriting £ in terms of £ using (j33|) and the expressions (|37h we have 

From this we see that, at long wavelengths, the leading solution £ oc 1/t corresponds to £ = const — 
the attractor we found studying homogeneous perturbations. Let us check that the same holds for 
the metric as well. In the absence of anisotropic stress ^ equals <3?|j. For the £ ~ 1/t mode, we can 
neglect the first term of the right-hand side of eq. (f53"j) so that 

* ~ -H£ with £ = const. (58) 

As £ = const describes the unperturbed solution shifted in time by an amount £, we see that the 
metric also describes the unperturbed FRW shifted in time by £. Indeed a 2 (t + £) — (! + 2H£)a 2 



What is quite confusing is that we are used in standard inflation to identifying the existence 
of an attractor with the £ = const mode's being the leading one at late times. Why does not 



°The mixing with gravity is important at all scales, but it fades away in the Mp\ —> oo limit with 
fixed / and Hq. Notice however that if one uses (|55|) to write the action for ip in the spatially flat gauge, 
£ = — ip ■ H/ttq, the result does not reduce to the one without gravity, eq. (|2Tj) . sending Mpi — > oo. The 
reason for this unexpected result is that if we take Mpi — > oo in a generic gauge, the spacetime does become 
flat, but we are not guaranteed that it be written in standard coordinates with metric 7/p„. Indeed it is 
straightforward to check, starting from (|32|) and (|33j) and doing the change of gauge, that in spatially flat 
gauge 

SN £-(£\ 2 C I 9MMt2 1 C (55) 

The limit must be taken while keeping (p, i.e. C/H , constant; we see that both SN and remain finite in 
spatially flat gauge when Afpi — > oo: the metric does not become rj^ u . This does not happen if one considers 
a model with a scalar field with a minimal kinetic term: in this case 5N and go to zero when Afpi — > oo. 
6 Actually this equality is not so straightforward to obtain. If one expresses $ as a function of £ using 
flSU) and (jnil), one gets 

•-4< + (££)'■ <m 

There is a partial cancellation between the two terms using the equation of motion of £ and the term H / H 2 Q 
cancels (notice that H /H 2 3> 1 in the limit we are considering.) To retain terms 0(1) x £ — which are needed 
to check whether $ = ^ — one should go beyond the approximation we are using. We do not do so, but we 
use the expression of $ in which no cancellation occurs, taking for granted that subleading corrections will 
enforce $ = *£. 

7 In order to cast the metric in FRW form, one should also do a redefinition of the time coordinate. It is 
easy to check that the resulting effect on the metric is suppressed by O(Ht) with respect to <3>. 
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this happen here? To understand this different behaviour, it is useful to follow [7] and study the 
most generic adiabatic mode in Newtonian gauge. One can do so by considering the residual gauge 
freedom of Newtonian gauge at k = : 

t^t + e(t) x i -W(l-A) , (59) 

with constant A. Under these transformations the potentials transform as 

9 -f^ + He- A $ ^ $ - e . (60) 

These are just gauge modes. To become the k — > limit of physical solutions, they must satisfy 
Einstein equations at infinitesimal but non- vanishing k. These will set $ = *$>, assuming we have 
no anisotropic stress. Using this, one can find the most generic adiabatic mode solving for e(t) 

e(t) = -A- f a(t')dt' + 4t . (61) 

V ' a{t) J V ; a{t) V ' 

A generic adiabatic mode is thus fixed by two constants: A and c. Usually one neglects c as e is 
dominated by the integral, which grows as t for any power-law expansion of the form a(t) oc t 2 ( 1 +' u )/ 3 
with 1 + w > 0. This is the standard adiabatic mode. Notice that the constant A is the value of £ 
for the constant solution [7J: indeed it parameterizes a rescaling of the spatial coordinates. For a 
constant w, one obtains a constant value for $ 

* = * = -l^X , (62) 

which is the standard relation between the Newtonian potential and the conserved quantity £. 
Our situation, however, is quite different. Taking a ~ 1, we see that the first term in eq. (|6ip 
still goes as t, but now this implies that it becomes subdominant as t gets close to zero. The 
dominant adiabatic mode is now given by e = c/a, which in our approximation amounts to e ~ c, 
corresponding to a constant time-shift. This makes perfect sense: as the model we are studying 
has a smooth Mp\ — > oo limit, we expect it to be dominated by an adiabatic mode that reduces to 
a constant time shift, which is what we get in the absence of gravity. 

It is immediate to work out the precise relation between the coefficients A and c, describing the 
general long wavelength solution in Newtonian gauge and the two modes of £, eq. (I49p . We can 
relate \& with £ using ([53|) and ([56]) : we see, as already discussed, that the ( oc cos kt/t corresponds 
to a constant e and has A = 0. On the other hand the mode £ = A sin kt/kt gives e ~ At with c = 
and A = A: as expected A corresponds to the amplitude of the constant £ mode. 

This adiabatic mode enjoys similar properties as the standard £ = const one. Independently of 
any details about the future evolution, e oc 1/a remains a solution. This implies that this mode is 
completely irrelevant for observations, as it quickly decays away in the standard FRW evolution. 

To complete the picture, in Appendix [A] we reproduce in Newtonian gauge the results for 
homogeneous perturbations of Section [2l 

5 A second scalar in the fake de Sitter 

The conclusion of the previous Sections is that the fluctuations of the ir scalar are cosmologically 
irrelevant. We are therefore forced to look for an alternative mechanism to give rise to the observed 
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scale invariant spectrum of primordial perturbations. We do not have to try very hard as the model 
itself naturally suggests one. The fictitious metric 

9® = e 2n{x) V,u , (63) 

with tt following the unperturbed solution (J2J), describes de Sitter space. Notice that any coupling 
of additional degrees of freedom with tt will have to go through the metric above to preserve the 
conformal invariance of the theory and, for de Sitter, any tensor constructed with the metric is 
proportional to the metric itself. This means that a second scalar a coupled to tt will behave as in 
de Sitter space. If a is massless — which can be ensured by a shift symmetry a — >• a + const — its 
spectrum will be scale invariant, 

(a(k)a(k')) = (2n) 3 6(k + k')^, (64) 

while the inclusion of a small mass term will tilt the spectrum either way. It is straightforward 
to check that the corrections coming from the evolution of the "real" metric g^ u are exponentially 
small, for modes of cosmological interest. This means that a massless a will acquire an exactly scale 
invariant spectrum, which is still marginally compatible with the data [ID]. We stress that gravity 
breaks the conformal symmetry of the Galileon Lagrangian, and this may induce small corrections 
to the scale invariant spectrum above. 

The conversion of a fluctuations into adiabatic ones can happen through one of the mechanisms 
that have been studied at length for inflation: a may change the way tt reheats the Universe 
[11[I12], or become relevant at a later epoch [13] . As the conversion mechanism is model dependent, 
unfortunately we cannot infer the value of Hq from the data. However, the experimental signatures 
are the typical ones for a "second field" mechanism: large, local non-Gaussianities and, possibly, 
isocurvature perturbations. 

As it is common to most alternatives to inflation, the gravitational wave spectrum is very blue 
in our model, and unobservable. Gravity waves are just sensitive to what the "real" metric is doing; 
given the rapid increase of H we will have a very blue spectrum of tensor modes. Indeed given that 
a can be approximated as a constant while H blows up, each mode gets frozen — k/a ~ H — at an 
amplitude of order of the Minkowski quantum uncertainty. Therefore we have a spectrum ex k^ 1 , 
very suppressed on cosmological scales. 

6 Faster than light 

We finally analyze the worrisome feature of our scenario: superluminality. It was shown in |15tl2|l3] 
that for DGP- and Galileon-like theories superluminal excitations may be generically expected 
about non-trivial solutions. Let's briefly review the general arguments of [2] and check whether 
they apply to our case as well. 

Consider our model (pi in the absence of gravity (Mpi — > oo, g^ v — > r/^ u ) and about the trivial 
configuration with tt = (jj). tt excitations are of course exactly luminal, for Lorentz invariance is 
unbroken. Now turn on localized sources so as to create a weak stationary ttq(x) field. By 'weak' 
we mean that 7r's self-interactions are unimportant to determine the field configuration, that is the 
solution obeys 

V 2 vr ~ , (65) 



8 As we discussed, fluctuations around tt = are ghost-like for the action (fT]), but this sign can be flipped 
starting from a more general Galileon Lagrangian and preserving the NEC violating solution we are interested 
in. Let us assume thus that the sign is the healthy one: the discussion below remains unaltered. 
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outside the sources. The quadratic Lagrangian for small perturbations 5tt about this solution is 

S 2 C = f G^ d^d7r d u 5ir + ..., (66) 

where the (inverse) effective metric is 

ff " = ^( 1 +5^ a »)-i^* ,8 ""»- (67> 

and the dots stand for non-derivative terms for 5tt, as well as for corrections that are at least 
quadratic in the background field 7To and derivatives thereof. The causal structure is determined 
by the highest derivative terms in the quadratic Lagrangian — therefore non-derivative ones like 
mass terms are irrelevant for our discussion. Also irrelevant is the correction proportional to rjau in 
(|67p — at the order in ttq we are keeping, it is just an overall conformal factor, which does not affect 
the light-cone aperture. In conclusion, the propagation of 5ir is constrained by the light-cone of 

4 1 
G^y ~ rj^ + --pjd^Avro . (68) 

Because of (|65p . this is narrower than the Minkowski light-cone in some directions, but wider in 
others [2] . Notice that this conclusion only relies on the presence of the Galileon cubic interaction 
(37r) 2 D7r, regardless of its sign — which in fact can be changed by redefining -k — > —it. 

At the classical level there is no way out — generic weak-field solutions admit superluminal 
excitations. At the quantum one however, we have to make sure that the effect be measurable 
within the effective theory. Roughly speaking, this is the case if signals of frequency lower than 
the UV cutoff can gain an order one phase-shift over exactly luminal signals of the same frequency. 
Suppose ddiTQ in (J68j) can be approximated as constant over some distance L, and let's make a 5tt 
signal travel such a distance. In order for our weak-field approximation to be valid throughout this 
region we need 

ddiT <C Hi , dvr (L) ~ ddn x L <C H , n (L) ~ ddn xi 2 «l. (69) 

The superluminal shift in the velocity is of order 5c ~ ddiro/HQ, corresponding to an overall 
phase-shift for 5-it 

<5phase ~ ScuiL ~ — — — , (70) 

Hq Hq 

where u is the signal's frequency. The phase shift is much smaller than oo/Hq in the weak field 
regime — see eq. (|69p — and becomes of order lo/Hq if we stretch the linear approximation to the limit. 
This means that the source of superluminality we identified (there may be others) is ineffective if 
we declare that the effective theory breaks down at energies/frequencies of order Hq, so that the 
would-be superluminal effect is not measurable consistently within the EFT. Notice that Ho is well 
below the strong-coupling scale of the theory [2]. 

The above discussion applies to small deformations of the trivial n = background. But what 
about our cosmological solution? As we will see, essentially the same conclusion applies. However 
the scale Hq will be replaced by 1/t — not surprisingly given our solution's scale-invariance. Consider 
first our solution in the absence of gravity, eq. ([2]). Given the high-degree of symmetry of the 
original Lagrangian (SO (4,2)) as well as of the solution (SO(4,l)), small fluctuations about this 
configuration are exactly luminal [2J . We now want to run the same argument as above — turning on 
a weak-field deformation of this background and studying its small excitations — but for simplicity 
we want to do so in a patch small enough so that we can approximate the "fake" de Sitter metric 
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g27r d s^ as £ a ^_ Given the homogeneity of de Sitter space, all points are equivalent, and we can 
choose for instance x = 0, t = —Hq as the center of our patch. We can then perform a special 
conformal transformation combined with a translation — recall that our Lagrangian is conformally 
invariant — to rewrite the de Sitter background as [TBI [3J 



1 
72 



e 2n *s = * , (71) 

l + ±H 2 x»x» 



where now x^ is measured from our new origin q This rewriting of the solution makes it immediate 
to carry out our analysis in a small patch centered at the origin. First of all, the de Sitter conformal 
factor reduces to 

e 2n ds „ 1 + 2^ „ 1 _ lfl-2 x ^ (72) 

Second, at lowest order in HqX 2 , the full non-linear dynamics of perturbations about the ir^s 
background are given by a Galileon Lagrangian whose coefficients are of the same order as the 
original ones [3]. In other words, expanding a generic Galileon Lagrangian about a ir oc x^x^ 
solution yields another Galileon Lagrangian with similar coefficients. As a consequence, as long as 
we restrict to distances from the origin smaller than Hq , our analysis above applies unaltered. 
We thus have superluminal excitations, which do not really have a chance of yielding measurable 
superluminal effects if the effective theory breaks down at frequencies of order Hq, or below. 

However here the appearance of Hq stems uniquely from our choosing to expand about t = Hq . 
This is made manifest by working with the field 

<P = H Q e* , (73) 

in terms of which the Lagrangian takes the form [2] 

£= i^*V'^rJ' (74) 

where F is a polynomial with order-one coefficients. The overall dimensionless factor f 2 /Hq has no 
effect at the level of classical equations of motion. Then the only scale present in the Lagrangian 
is the local value of <ft. Our de Sitter solution (in the original coordinates) corresponds to 

0ds = ~\ ■ (75) 

So, the fact that to avoid superluminality in a neighborhood of t = Hq we have to impose a 
frequency cutoff of order Hq, implies that to avoid superluminality about a generic t the UV cutoff 
has to be l/|i|. Notice that cosmological perturbations of the Galileon field freeze-out at frequencies 
precisely of order l/\t\. Therefore, if we decide to ban superluminality from our effective theory, we 
also lose predictivity for cosmological observables. The other possibility — swallowing the presence 
of measurable superluminal effects in a Lorentz-invariant effective theory — does not necessarily 
lead to inconsistencies. As long as the effective theory is free of closed time- like curves (which for 
our model has not been proven nor disproven yet), there are no pathologies from the low-energy 
viewpoint (see for example |14] for an optimistic point of view). However, physically measurable 
superluminality certainly implies that the effective theory at hand cannot arise as the low-energy 



9 Notice that performing conformal transformations does not perturb the causal structure, even though 
they do not commute with the Lorentz group, because they only affect the metric by a conformal factor. 
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limit of a microscopic theory with the standard relativistic causal structure, such as a renormalizable 
Lorentz-invariant QFT for instance |15j . 

The introduction of dynamical gravity perturbs the above analysis, and to some extent its 
conclusions too. To begin with, our cosmological solution ([2]) gets modified: slightly at early times 
(eq. (fT3j) ). drastically at late ones (eq. (HU)). As a consequence its de Sitter symmetry is gone, 
and small 5ir perturbations are no more exactly luminal, even in the absence of the sources we 
needed to run the above arguments. Second, given the peculiar structure of the Galilean self- 
interactions, the mixing of 5ir with scalar gravitational perturbations is relevant at all scales (see 
a related discussion in [5]), and cannot be ignored even when studying sub-horizon perturbations. 
We thus have to use eq. (|34h . which is the quadratic action for the propagating scalar mode about 
the FRW background, taking into account gravitational corrections. The propagation speed for 
short-wavelength excitations as measured by a comoving observer in terms of the background FRW 

metric is 

9 Bit) , N 

C C = TTT • ( 76 ) 



In other words, if c? thus defined is larger than one (, excitations exit the FRW light-cone. Plugging 
the approximate solutions we found in sect. [2] into eqs. ([351136]) . we get 

32 f 2 1 

^-°° : t^-gT&W (77) 

t ->> t : 4 ~* ° • ( 78 ) 

In both regimes the correction to the propagation speed is suMuminal — extremely so at late times 
l 10 l . This relaxes our conclusions above somewhat: since the cosmological background introduces a 
subluminal offset into the excitations' speed, we now need perhaps small, but finite perturbations 
to overturn this offset and make excitations about a new background super luminal. Although this 
is certainly more welcome than the opposite result, in practice it is not very helpful: at very early 
times the gravitational correction to the propagation speed goes to zero — like all other gravitational 
effects in our model. This means that, at least at early times, we have to live with superluminality, 
or give up the model. 

7 Conclusions and Outlook 

We are putting forward a model for our Universe's early cosmology that departs strikingly from 
the conventional inflationary picture. Schematically: there is no Big Bang in our past; spacetime 
is flat at t — ¥ — oo; related to this, the Universe is initially devoid of any form of energy; energy 
and the associated Hubble expansion get created by a NEC-violating sector. 

The main virtue of the model lies not in its sheer radicalness — we are certainly not the first 
authors to come up with "phantom" -like equations of state — but in its being able to associate 
such a radicalness with an healthy effective field theory coupled to gravity, and in the consequent 
robustness of the scenario. Our theory — the Galileon or more precisely its conformally invariant 
generalization — is well-behaved classically as well as quantum mechanically, even for strongly non- 
linear background solutions. The structure of the Lagrangian is protected by symmetries. More 
relevant for us, the system retains stability even when the stress-energy tensor violates the NEC 



10 The vanishing of c? for t approaching to signals that, in such a limit, higher-derivative corrections to 
the perturbations' gradient energy cannot (and should not) be neglected, pretty much like for the ghost 
condensate [4]. 
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[2]. Partially as a consequence of this, the cosmological solution we outlined is an attractor — the 
universe wants to follow it even when initially displaced from it — which implies that our model 
solves the horizon and flatness problems as well as standard inflation does. Remarkably, expansion 
is not put in as an initial condition but follows from generic initial conditions, including contracting 
ones — within a bounded basin of attraction of course. As to density perturbations, in its minimal 
incarnation our model does not produce sizable adiabatic perturbations on cosmological scales. 
However the symmetry structure is so constraining that postulating the existence of extra light 
scalars unavoidably yields nearly scale-invariant spectra for them — which can later be converted into 
adiabatic perturbations via any of the standard conversion mechanisms available on the market. 
The downside is that, like for standard multi-field inflationary models, predictions are more model- 
dependent than for single-field slow-roll inflation, although the presence of a sizeable local Non- 
Gaussianity is rather robust. 

The only reservation we have about welcoming our model as a compelling alternative to infla- 
tion concerns superluminality. The generic presence of superluminal excitations about non-trivial 
solutions indicates that our model cannot arise as the low-energy limit of a standard relativistic 
UV-complete theory, like e.g. a renormalizable Lorentz-invariant QFT. Depending on one's personal 
taste and attitude, reactions to giving up Lorentz invariance in the UV may range from disgust to 
excitement. Minimally, it is fair to say that it makes us depart from known territory, especially 
when gravity is involved. We therefore feel that it deserves special care. It is interesting to note 
that superluminality in our model is tied to the presence of the DGP-like interaction ((97r) 2 D7r, 
which in turn is forced upon us by demanding that scattering amplitudes obey standard properties 
of S'-matrix theory [2]. However it is possible, in principle, that suitable deformations of the theory 
exist where such an interaction is absent and where superluminality is gone as wellllj Of course we 
would like to maintain the nice features of our cosmological scenario — among which the near de- 
Sitter invariance of the solution. So, instead of considering a conformal completion of the Galileon 
theory, we may consider a different symmetry group containing the de Sitter one as a subgroup and 
that reduces to the galileon symmetry group in the appropriate limit. An obvious choice is the 5D 
Poincare group ISO(4,l) [31 El [20]. This possibility certainly deserves further study. 
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A Details on adiabatic perturbations 

This Appendix complements the study of adiabatic perturbations of Sections [3] and HI 



11 Given the non-renormalization theorem of [19] . such a "tuning" would be preserved by quantum correc- 
tions. 
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Unitary gauge action. To deduce eq. (J26J), we use that in unitary gauge the various terms in 
the action read 



-ge 2 *(dvf 

^nvr(a^) 2 



■ 4 ,- 1 



1 
N 



N 3 

-27T^7f Vff3^3 + *0 V93JJ 



N l di 



N 2 



dt 



N 2 



(79) 
(80) 
(81) 



where g% = det gij, we used that \J—g = N^fg^, and in the last line we integrated by parts. Also 
we made use of g° l = N l /N 2 , where and henceforth spatial indices are raised and lowered with 
the spatial metric gij . We can rewrite the terms in bracket in terms of the extrinsic curvature of 
costant-t hyper surf aces, 

K i3 = -^ [d tgij - ViNj - V 3 Ni] . (82) 

Indeed after straightforward manipulations and integrating by parts we get 



^0 V93j-f 



N*di^ - a-* 



] N 2 



N 2 



2*0*0 V53T73 + 3*0 V93TH K 



'AT3 



N 2 ' 



(83) 



where we used that dt ^/gi = n^/s E g 1 ^ dtgij • The first piece cancels exactly the first term in eq. (|8ip 
and we are left with eq. (|26p (j I ). 

Unitary action in the standard form. To cast eq. (|26p in the form (|28p we can expand the third 
term of (|26p as 1/iV 4 = 2/N 2 — 1 + 4<5 N 2 + . . . The second term can be rewritten as 



±K! = 6±5Ki + Ki + 3H±-3H 



(84) 



Notice that K\ appears as an additional tadpole term besides the ones in eq. (J28J). However one 
can get rid of it using the identity [17J 13 I 



<?xV=gf(t)K£ = / d*xV=gf(t)V^ 



d 4 xV^5 M /(t)^ 



d^Xy 



-9f(t)^. 



(85) 
In this way one can write the action in the form (|28p and check that the coefficients of the tadpole 
terms can indeed be written in terms of H and H using the expression of the ir stress-energy tensor. 
Constraint equations. The explicit form of the constraint equations is 



M 2 ! 



Ri 



' 'EijE'i -E\ 2 ) + ^H-2(3H 2 + H) 



N 2 



N 2 ' 



+ 2M*8N - 2M 6 8E % 



V,. 



1 



M 2 !- (E l j - 5) E k k ) - 5) MHN 



l N 











(86) 

(87) 



whose solution at linear order gives (I32p and (I33p . 

Homogeneous perturbations in Newtonian gauge. Let us check that in Newtonian gauge we find 
the same two regimes of perturbations we found in Section [2) one when the perturbation dominates 



12 As tttt 2 is a total derivative, each term in (<97r) 2 D7r contains at least two spatial derivatives. That's why 
it is not surprising that it can be written solely in terms of an operator containing the extrinsic curvature 
K, which contains two spatial derivatives on n. 

13 Notice there is a sign error in the last term of eq. (80) of [T7] . 
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the energy density (this is equivalent to decoupling gravity as the energy density of the background 
vanishes when Mp\ — > oo) and one when perturbations are small and only perturb the background 
energy density. The first regime is obtained simply by sending Mp\ — > oo. In eq. ()56p the leading 
solution in £ gives, as we saw, the £ = const mode, while the solution £ oc sinkt/kt gives the 
decaying solution £ ~ t 5 , as it can be easily verified. Actually this identification holds not only on 
large wavelengths, but it is exact in the limit Mp\ —> oo: the second mode of £ in eq. (J49H matches 
with the solution of eq. (|2ip that is constant at small k, while the first one matches with the £ mode 
that as t 5 at small k, without mixing. It is much trickier to study the regime when perturbations 
do not dominate the energy density. From Section [2] we expect the decaying mode to give £ ~ t 2 
in this case, but as we said the decaying mode of C gives £ ~ t 5 . The trick is that one has to be 
careful about the two limits Mp\ — > oo and k — > 0. Indeed, if one keeps higher orders in 1/Mpj in 
the action for (, eq. (J39j) . in the relation between £ and £, eq. (|56|) . and takes also into account the 
change in the time variable to compare with the FRW solutions of Section [2j the decaying mode 
goes as 

+ H 2 t 5 k 2 o / 13vr 13 2l \ 

^^7^ M Pl+ -Tr + Tir'rt* +■■■ ( 88 ) 



5/ 2 ^ V 5Jfe 15 

where the dots stand for terms of higher order in 1/Mp, and k. We see that the £ ~ £ 5 solution 
dominates if one sends Mp\ — > oo at fixed k. On the other hand sending k — > at fixed Mp\ we 
have a constant term, which describes a mixing with the dominant £ = const mode and the t 2 
behaviour we were looking for. 

B Squeezing and absence thereof 

As we discussed, the perturbations of the Galileon field tr have quite peculiar properties, rather 
different from the standard inflationary scenario. The best way to pin down the status of these 
perturbations is to determine the quantum state of each Fourier mode, when it comes back into the 
Hubble radius, for example during a radiation dominated phase. Indeed, in the linear approxima- 
tion, the field is just a collection of harmonic oscillators, each with a time dependent Lagrangian, as 
the background we are perturbing around is time dependent. This causes each Fourier mode to be 
in a squeezed state, when it gets back into the Hubble radius. The amount and direction of squeez- 
ing uniquely fix the state of the perturbation and also tell us whether a classical interpretation in 
terms of classical stochastic variables is possible. 

In order to follow the evolution of each Fourier mode until the it comes back into the Hubble 
radius, it is useful to have a unified description in which perturbations around a homogeneous, 
isotropic Universe are always described in each phase of the evolution of the Universe by an action 
for the same scalar variable. We will do so using £ as such a variable. In this Appendix we start 
calculating the squeezing status of perturbations in the case of inflation and then compare it to our 
Galileon case. Of course the case of inflation is quite well known, but the way it is presented here 
is, to our knowledge, new. 

The action for £ during inflation is of the form [6] 



'Pi 



S = Ml / d A x a 3 e 



C 2 -^(VC) 2 



(89) 



where e = <p 2 / \2H 2 -M™) can be taken as a constant at leading order in slow-roll. The field C, can 
be decomposed in terms of annihilation and creation operators as 
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As the field satisfies the equations of motion, the functions Cr (t) are solutions of the equations of 
motion that at very early times, when the mode is much shorter than the Hubble radius, reduce to 
the Minkowski form: 

Cf(t) = -=L %= (l - z-M ^ +fkS • (91) 

Let us define a scalar product between classical solutions of the field equation 

(Cf; c 2 cl ) = -if d 3 x (cfnf - cfn? ) , (92) 

where II is the momentum conjugate to £. For the action (189h . II = 2Mpia 3 e£ . It is important 
to notice that this scalar product is time independent as a consequence of the equations of motion. 
The solutions £ are normalized as 

(Ct;Ct) = -(^f6(k-k f ) (93) 



5 fc ' *k' 

ACl. ACl 



(c?-x?;) = o. 



The action describing scalar perturbations during a phase dominated by a barotropic fluid with 
p = wp is given by [18] 



> .l/iW^« 3 ^ 



a z 



C 2 -3( V 2 • (94) 



Let us concentrate for example on a period of radiation dominance, w = 1/3, which gives an 
evolution a oc t 1 ' 2 , H = l/(2£). One can still perform an expansion analogous to (|9Up . 

c(t ' f) = / Wf &m + cf (t)at) . (95) 

Now the appropriate solutions of the equation of motion are given by 

C ~£i( t ) = _5 J: . V^+^ . ( 96 ) 



12M PU/2fc/V3 a 



These functions reduce to the Minkowski result in the limit in which the modes are well within the 
Hubble radius and they are normalized in the same way as in (|93j) . The choice of phase of these 
solutions is done for later convenience. 

Equating the two expansions (1901) and ([95]) for each Fourier mode we get 



(far + C cl ra f - = Cfa t + C cl * r a f - (97) 

The two set of modes are related by 

Cf = « j $ + P U C% (98) 

where the coefficients can be generically calculated using the scalar product! 14 ! 

«*=<?;#> ^ = -(^0> (99) 



14 Notice that the scalar product (|92[) is well defined at any time, and time independent. Of course the 
explicit expression for the momentum II in terms of £ depends on the action for £ valid at any given moment. 
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but at small k we can use a quicker method — see below. Plugging eq. (j98|) into (J9TJ) gives the 
relation among the two sets of creation and annihilation operators 

H = ^H + ^U- ( 10 °) 

We assume to be in the vacuum at the beginning, i.e. in a state which is annihilated by ar. This 
state does not evolve in time as we are putting the time evolution in the operators. In terms of 
the "radiation dominance" operators aV, the state of each harmonic oscillator is annihilated by a 
linear combination of ar and a-: it is a squeezed state. Taking the norm of both sides of eq. (|98p 
we derive that the coefficients a and (3 must satisfy 

l«/ " l/%| 2 = 1 ( 101 ) 

which is equivalent to the condition that both sets of creation and annihilation operators satisfy 
the standard commutation rules. 

The simplest way to relate the modes £ and £ , is noticing that their real parts become 
constant in the long wavelength limit, while their imaginary parts have constant conjugate momenta 
in the same limit. Given the independence or incompatibility of these two conditions — constant 
field vs. constant momentum — this implies that real and imaginary parts do not mix in the long 
wavelength limit, that is 

CS 1 + Cf = A k (Cf + Cf ) , CI 1 - Cf = B k (Cf - Cf) , k^O. (102) 

Condition (jlOip implies that the two coefficients A k and B k are the inverse of each other: A k = B^ . 
A k and B k can be calculated using the explicit expression of the modes: it is enough in the two 
cases to compare the value of C, and II which is approached in the long wavelength limit. One gets 

A k = #l*Jl «« H « H ™ Bk = l_. (103) 

These expressions are time-independent, as they should. It is easy to realize that, neglecting 
numerical factors, one has 

A k ~ — > 1 , (104) 

flout 

the ratio between the scale factors when the modes leave the Hubble radius and when they come 
back in. The uncertainty in the real part of (, is huge compared to that we would have in the 
"radiation dominance vacuum" . Conversely the uncertainty in the imaginary part is very suppressed 
with respect to the vacuum state. If we neglect this minuscule uncertainty in the imaginary part, 
assuming that fluctuations in every observable will be dominated by the huge real part's uncertainty, 
we can treat the quantum field as a classical stochastic variable. 

Let us now come back to our model, and carry out exactly the same procedure. Now the £2 
modes must be calculated from the action (|39p and they are given by 

(£= * ]= l e -ikt+ika! ^ (105) 



'/>• 



3y/2M^H V2k t 



while the radiation dominance modes £r remain, obviously, the same. Again the real part of the 
modes gives £ = const in the long wavelength limit, while the imaginary part has II = const. 
Matching these constants allows to calculate the coefficients A and B 

A - - Mhr^ H S B ^i b - <m > 
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Let us assume that the Galileon dominated phase ends when H ~ Mpy/fHQ, i.e. when the approx- 
imation of treating gravity as a perturbation breaks down, and that immediately after we have a 
radiation dominated phase. Evaluating A k at the transition between the two phases we get 

A k ~ B k ~ 1 . (107) 

There is no relevant squeezing of the modes, so that during radiation dominance the uncertainties 
are close to the standard zero point quantum fluctuations. Of course in this case the perturbations 
are completely negligible and no classical interpretation is possible 15 !. 

A sizable generation of perturbations needs a large squeezing. This, in inflation, is closely 
related to the existence of a dynamical attractor which makes the C, = const mode dominate. Our 
model shows that the existence of an attractor does not by itself guarantee a sizable generation of 
perturbations: during radiation dominance the modes are essentially in their vacuum state. 
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